A new method to measure the attenuation of hadrons in extensive air showers 
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Extensive air showers are generated through interactions of high-energy cosmic rays impinging 
the Earth's atmosphere. A new method is described to infer the attenuation of hadrons in air 
showers. The numbers of electrons and muons, registered with the scintillator array of the KAS- 
CADE experiment are used to estimate the energy of the shower inducing primary particle. A large 
hadron calorimeter is used to measure the hadronic energy reaching observation level. The ratio of 
energy reaching ground level to the energy of the primary particle is used to derive an attenuation 
length of hadrons in air showers. In the energy range from 10^ GeV to 3 ■ lO'^ GeV the attenuation 
length obtained increases from 170 g/cm^ to 210 g/cw? . The experimental results are compared to 
predictions of simulations based on contemporary high energy interaction models. 
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ergies achieved in man made accelerators by orders of 
magnitude. Hence, in the literature many attempts are 
described to extract properties of high-energy hadronic 
interactions from air showers induced by cosmic rays in 
the atmosphere. Among the most interesting quantities 
is the attenuation length of hadrons (e.g. Q p. 162), in 
theoretical considerations closely connected to the inelas- 
tic cross section. 



In the present work, we use the energy absorbed in 
a material within a certain atmospheric depth X to de- 
fine an attenuation length. In this new approach we use 
the number of electrons and muons, registered with a 
detector array to estimate the energy of the shower in- 
ducing primary particle Eq, see (|4]) below. The energy 
reaching the observation level in form of hadrons ^ Eh 
is measured with a hadron calorimeter. The fraction of 



I. INTRODUCTION 



Since the earliest days of cosmic-ray investigations it 
has been realized that these particles provide a unique 
possibility to study interactions at high energies [l|, 0. 
Even today the energies of cosmic rays exceed the en- 
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surviving energy in form of hadrons is defined as 



R = 



En 



The attenuation length A^; is then defined as 



or 



TiEh — Eq exp 



R = exp — - — 



X 



(1) 



(2) 



(3) 



In the literature different definitions for the attenua- 
tion length in air showers are introduced. Frequently, the 
attenuation length is derived from measurements of the 
electromagnetic shower component, e.g. [1, 0, H, 0, H, 0, 
[lOj . Investigating single hadrons the attenuation length 
is related to the absorption of hadrons in the atmosphere 
Pioneering work to derive inelastic cross sections 
from the measurement of single hadrons has been con- 
ducted by Yodh and colleagues , later followed by 
the prototype of the KASCADE calorimeter [l3[. 

The new approach presented here is complementary to 
the different methods described in the literature. In con- 
trast to methods using the electromagnetic shower com- 
ponent, the present work focuses directly on measure- 
ments of hadrons to derive an attenuation length for this 
shower component. The values obtained are not a pri- 
ori comparable to other attenuation lengths since they 
are based on different definitions. It should be noted 
that the experimentally obtained attenuation length is 
affected by statistical fluctuations during the develop- 
ment of the showers. However, in the present work we 
do not attempt to correct for this effect. 

After a description of the experimental situation 
(Sect.|lT|), the experimental results and comparisons with 
air shower simulations are described in Sect.lIIll 



II. EXPERIMENTAL SET-UP 

A. The apparatus 

The experiment KASCADE, located on site of the 
Forschungszentrum Karlsruhe, 110 m a.s.l., consists of 
several detector systems. A description of the perfor- 
mance of the experiment can be found elsewhere [14| . A 
200 X 200 m^ array of 252 detector stations, equipped 
with scintillation counters, measures the electromagnetic 
and, below a lead/iron shielding, the muonic components 
of air showers. In its center, an iron sampling calorime- 
ter of 16 X 20 m^ area detects hadronic particles. The 
calorimeter is equipped with 11 000 warm- liquid ioniza- 
tion chambers arranged in nine layers. Due to its fine 
segmentation (25 x 25 cm^), energy, position, and an- 
gle of incidence can be measured for individual hadrons. 



A detailed description of the calorimeter and its perfor- 
mance can be found in [TH , it has been calibrated with a 
test beam at the SPS at CERN up to 350 GeV particle 
energy jl6| . 



B. Observables and event selection 

The position of the shower axis and the angle of in- 
cidence of a cascade are reconstructed by the array de- 
tectors. The total numbers of electrons iVe and muons 
are determined by integrating their lateral distribu- 
tions. In case of muons, the truncated muon number 
N'^ is used for experimental reasons. It is the number of 
muons integrated in the distance range 40 — 200 m from 
the shower axis. For a detailed description of the recon- 
struction algorithms see [iTj . The position of the shower 
axis is reconstructed with an accuracy better than 2 m 
and the angle of incidence better than 0.5°. 

The hadrons in the calorimeter are reconstructed by 
a pattern recognition algorithm, optimized to recognize 
as many hadrons in a shower core as possible. Details 
can be found in [l3]- Hadrons of equal energy can still 
be separated with a probability of 50% at a distance of 
40 cm. The reconstruction efficiency rises from 70% at 
50 GeV to nearly 100% at 100 GeV. The energy resolu- 
tion improves from 30% at 50 GeV to 15% at 10"^ GeV. 
The hadron number Nh and hadronic energy sum ^ Eh 
are determined by the sum over all hadrons in a distance 
up to 10 m from the shower axis. A correction for the 
missing area beyond the boundaries of the calorimeter is 
applied. The hadron lateral distributions are relatively 
steep, the hadronic energy density decreases by about 
two orders of magnitude within the first 10 m from the 
shower axis [l^. Therefore, it is sufficient to measure 
hadrons in a relatively narrow range around the shower 
axis only in order to collect a significant fraction of the 
total hadron energy. The observable ^ Eh includes also 
energy of hadrons which could not be reconstructed in- 
dependently, because they are too close to each other. It 
shows up in the simulated and experimental data in the 
same manner. 

To be accepted for the present analysis, an air shower 
has to fulfill several requirements: at least one hadron 
has been reconstructed in the calorimeter with an energy 
larger than 50 GeV, the shower axis is located inside 
the calorimeter, the electromagnetic shower size Ne is 
larger than 10^, the truncated muon number N'^ is larger 
than 10^, i.e. the primary energy is greater than about 
3 • 10^ GeV, and the reconstructed zenith angle is smaller 
than 30°. From May 1998 to October 2005 312000 show- 
ers have been measured meeting the criteria mentioned. 

To avoid corrections for different angles of incidence 
the following analysis is restricted to showers with zenith 
angles 9 < 18°. The primary energy Eq of the shower 
inducing particle is roughly estimated based on the num- 
ber of electrons and muons N'^^^ registered with the 
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FIG. 1: Fraction of energy TiEh/Eo reaching the ground in 
form of hadrons as function of estimated primary energy Eo 
for all data and for selections of light and heavy primary par- 
ticles. 
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KASCADE scintillator array 



\gEo « 0.191giVe + 0.79 IgiV;, + 2.33. 



(4) 



The average ground pressure during the observation time 
amounts to 1004 hPa, corresponding to an average at- 
mospheric column density Xq—1Q2'6 g/cm^. The attenu- 
ation is measured at this depth Xq, the average vertical 
thickness of the atmosphere above the KASCADE exper- 
iment. 



C. Simulations 

The shower simulations were performed using COR- 
SIKA [l^. Hadronic interactions at low energies were 
modeled using the FLUKA code [13, HI. High-energy 
interactions were treated with QGSJETOl [22| {E > 
200 GeV). Showers initiated by primary protons as well 
as helium, carbon, silicon, and iron nuclei have been 
simulated. The simulations covered the energy range 
10^ - 10® GeV with zenith angles in the interval 0° - 32°. 
The energy distribution of the showers followed a power 
law with a spectral index of —2.0. For the analysis the 
energy distribution was converted to a power law with 
an index of —2.7 below and —3.1 above the knee with 
a rigidity dependent knee position (3 • 10^ GeV for pro- 
tons). The positions of the shower axes are distributed 
uniformly over an area exceeding the calorimeter surface 
by 2 m on each side. In order to determine the signals in 
the individual detectors, all secondary particles at ground 
level are passed through a detector simulation program 
using the GEANT package [23]. In this way, the instru- 
mental response is taken into account and the simulated 
events are analyzed by the same code as the experimen- 



tal data, an important aspect to avoid biases by pattern 
recognition and reconstruction algorithms. 



III. RESULTS 

A. Surviving Hadronic Energy 

The energy of the primary particle is estimated from 
measurements of the number of electrons and muons in 
the shower with the scintillator array, see ([4]). The sur- 
viving energy in form of hadrons S-E/f is measured with 
the hadron calorimeter. A fraction R, see ([IJ of hadronic 
energy reaching ground level can be inferred as function 
of primary energy, as shown in Fig.[T] All error bars rep- 
resent statistical uncertainties only. Below 10^ GeV the 
values are affected by reconstruction efficiencies. In par- 
ticular, showers induced by heavy elements are less likely 
to be registered. Therefore, values are shown only for en- 
ergies exceeding 10^ GeV. Above lO"^ GeV the flux of the 
light cosmic-ray component decreases and the composi- 
tion becomes more and more heavy [23] . Mos likely, this 
causes the structures seen in the figure for energies ex- 
ceeding 10^ GeV. In the energy range investigated about 
0.3% to 0.8% of the primary energy reaches the observa- 
tion level in form of hadrons, most of them being pions 

In the energy range of interest the elasticity of pions 
depends only weakly on energy and can be approximated 
as e « 0.25 to 0.3 ]2^. With the relation R = e^, the 
average number of generations N in the shower can be 
estimated and it turns out that the registered hadrons 
(with energies above 50 GeV) have undergone about four 
to five interactions only. This number is confirmed by full 
air shower simulations. 

The fraction of hadronic energy reaching observation 
level increases with energy, since the effect of deeper pen- 
etrating showers clearly dominates over the small effect 
caused by the increase of the inelastic cross sections. 

The two-dimensional distribution of the number of 
electrons and muons for the measured showers is depicted 
in Fig. [21 The asterisks represent the most probable val- 
ues of the distribution. Also simulated N„ 



N' distri- 

butions have been investigated for primary protons, as 
well as helium, carbon, silicon, and iron-induced show- 
ers. Examples for protons and iron are depicted in Fig.[3l 
The solid lines represent fits to the most probable values 
represented by the circles and squares, respectively. It 
turned out that the slopes of the fits of all elements are 
about equal. The fitted line for carbon, parameterized as 



lg(7Ve) = 1.341g(iV')-0.15 



(5) 



is used in the following to divide the data into a sam- 
ple induced by "light" and "heavy" primary particles, 
respectively. It is indicated in Figs. [2] and [3] as dashed 
lines. In Fig.[3]it can be seen that (O indeed separates the 
data set into "light" and "heavy" . Almost no iron show- 
ers are above the dashed line and only a small fraction 
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FIG. 2: Number of electrons and muons for measured showers 
with zenith angle & < 18°. The most probable values of 
the distribution are indicated by the asterisks, the solid line 
represents a fit to this data. The dashed line represents ([5]). 

of proton induced showers is below the line. Especially 
the most probable values for protons and iron-induced 
showers are clearly above and below the dashed line, re- 
spectively. Parameterization ([5]) almost coincides with 
the most probable values of the measurements. 

Applying the selection criterion (O to the data, the en- 
ergy fraction reaching observation level is shown in Fig.[T] 
as well for light and heavy primaries. As expected from 
a simple superposition model, proton-like showers pen- 
etrate deeper into the atmosphere and transport more 
energy to the observation level as compared to iron-like 
showers. 



B. Attenuation Length 

Using the energy absorbed in the atmosphere at a 
depth Xq, the attenuation length Xe has been derived 
from the measured energy fraction, see ([3]). The results 
are presented as function of the estimated primary en- 
ergy in Fig. [5] a). Values for the complete data set as well 
as for the light and heavy selections are shown. The val- 
ues are compared to results obtained from full air shower 
simulations for primary protons and iron nuclei using the 
CORSIKA program with the hadronic interaction gener- 
ator QGSJETOl. The cuts for a "Hght" and "heavy" 
component were applied to the proton and iron simula- 
tions in the same way as for the measurements. 

The values shown in Fig. [5] have been obtained by ap- 
plying the same (quality) cuts and reconstruction algo- 
rithms to measured and simulated data. Thus, many 
uncertainties are expected to cancel. For the remaining 
differences between simulated and measured data a sys- 
tematic error for the hadronic energy sum Eh) and 
the total energy (Eq) of 10% each is assumed. With ([T]) 
this results in a 14% uncertainty of R. In turn, follow- 
ing ^ yields a systematic error for As of order of 2%. 



Therefore, the values shown in Fig. [4] have a systematic 
uncertainty of order of 2 to 4 g/cm^. As discussed above, 
the structures above an energy of 10^ GeV are most likely 
due to the chaning composition as function of energy. 

A comparison of the "heavy" selection with iron- 
induced showers shows that the "heavy" selection lies 
slightly above the simulated values. This makes sense, 
since the measured data contain a mixture of many el- 
ements, most of them being lighter than iron, thus the 
measurements should be above the iron points. On the 
other hand, looking at the "light" selection compared 
to the proton points from the simulations, one recognizes 
that at high energies data points are above the simulated 
values. This cannot be explained by a mixed composition 
and may be a hint towards a problem in the hadronic in- 
teraction model QGSJETOl. A possible explanation is 
that the attenuation length is too small, i.e. the cross 
section is too large. 

To test this hypothesis simulations have been car- 
ried out with a modified version of QGSJETOl, namely 
model 3a in Ref. [1^. The inelastic hadronic cross sec- 
tions have been lowered, e.g. the proton-air cross section 
at 10^ GeV is reduced by 5% from 385 mb to 364 mb 
and the elasticity has been increased by about 12%. A 
similar trend to lower cross sections has been found as 
well by the E AS-TOP experiment, with a value of 

at = 2 TeV (w 2 • 10^ GeV) At the high- 

est energies the lower proton-air cross section (443 mb 
at 10^ GeV) is compatible with recent results from the 
HiRes experiment 

<-Lr = 456 ± 17(3,^,) + 39(,y3) - ll(3y3) mb 

at 3 • 10^ GeV Q . The lower cross sections have been pro- 
posed originally to reduce the discrepancy in the mean 
logarithmic mass derived from experiments observing 
shower maximum and investigating particle distributions 
at ground level [2^, [2^ . Applying the altered version of 
QGSJET also slightly modifies the number of electrons 
and muons predicted at ground level. At energies around 
the knee (« 4-10^ GeV) the number of electrons increases 
by about 5% and the number of muons rises by about 

15% m. 

The corresponding results for Xe are presented in Fig.|^ 
b) . The typical energies of particles in a shower induced 
by a heavy nuclei are smaller than the typical energies 
in a proton-induced shower of the same primary energy 
(superposition model of showers). Thus, the effect of the 
modifications is stronger for proton-induced showers in 
Fig.dlb), see also [25j. In contrast to Fig. 2] a) the sim- 
ulated points for protons are now above the data points 
for the "light" selection. As a result, with the modified 
version of the interaction model an overall improvement 
of the situation has been achieved. 

It should be pointed out that the experimentally ac- 
cessible attenuation length A^; is extremely sensitive to 
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FIG. 3: Number of electrons vs. number of muons in simulated air showers for primary protons (a) and iron nuclei (b). The 
solid line indicates a fit to the most probable values (circles and squares, respectively), the dashed lines represent ([5]). 
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FIG. 4: Attenuation length Xe as function of estimated primary energy. The light and heavy groups in the measurements are 
compared to simulations for primary protons and iron-induced showers using CORSIKA with the hadronic interaction model 
QGSJETOl (a) and a modified version with lower cross sections and higher elasticity (b, model 3a in Ref. [25|). 



the inelastic hadronic cross sections. A relatively small 
modification {ap-air is changed by 5% only at 10^ GeV) 
yields significant changes, as can be inferred from Fig.lH 

In addition, the results have been calculated also as 
function of the hadronic energy sum at observation level. 
The results for all data, as well as for the light and heavy 
selections are presented in Fig. [5] a). Again, a closer in- 
spection yields unreasonable results when the "light" ele- 
ments are considered. At high energies even all measured 
events are above the proton simulations. In agreement 
with the previous discussion, the issue can be resolved by 
introducing a modified version of the interaction model, 
as can be seen in Fig.[S]b). Arranging the data in Y.Eh 
bins implies an enrichment of light primaries, which ex- 



plains why all data almost agree with pure proton simu- 
lations and the discrepancy between measurements and 
QGSJET 01 predictions are magnified. 



IV. SUMMARY AND CONCLUSIONS 

A new method has been developed to derive the at- 
tenuation length of hadrons from measurements of high- 
energy cosmic rays interacting in the Earth's atmosphere. 
The fraction of the energy of the primary particle reach- 
ing ground level in form of hadrons in air showers has 
been measured with the KASCADE experiment to in- 
crease from about 0.3% at 10^ GeV to 0.8% at S-IO"^ GeV. 
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FIG. 5: Attenuation length Xe as function of the measured hadronic energy sum at observation level. The light and heavy 
groups in the measurements are compared to simulations for primary protons and iron-induced showers using CORSIKA with 
the hadronic interaction model QGSJETOl (a) and a modified version with lower cross sections and higher elasticity (b, 
model 3a in Ref. |25ll). 



An attenuation length based on the absorbed energy has 
been defined. Corresponding values increase with energy 
from about 170 g/cm^ to « 210 g/cm^. 

A closer inspection of the attenuation lengths obtained 
for showers induced by "light" and "heavy" elements in- 
dicates that the cross sections in the hadronic interaction 
model QGSJETOl may be too large and the elasticity 
may be too small. A modification with altered parame- 
ters improves the situation. 

As final remark, it should be pointed out that the at- 
tenuation length As is extremely sensitive to inelastic 
hadronic cross sections. The sensitivity of air shower 
measurements, in particular of the hadronic component 
to properties of hadronic interactions has been demon- 
strated previously. For example, the dependence of 
observable quantities on the transverse momentum in 
hadronic interactions ^2T\ or on low-energy inelastic cross 



sections [28|. In a similar way, the data presented and 
the method introduced in the present article may serve to 
check and improve further hadronic interaction models. 
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